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PAaE  1«. 

The  far  right  of  Eq.  (33)  should  read  M  (H)  not 
M  (H).  ^ 

PAflE  19,  line  11, 

K  is  -  I30G  not  130a. 

H 

PAGE  19. 

Eq.  (37)  should  be  gpH  ^  =  gpH  -  4k  «  +  4tt  gPM  -  4k 

tot  ya  19  ^ 

= .309  cm.-^ 

Since  all  the  exchange  Integrals  are  negative,  all  numeri¬ 
cal  values  quoted  should  be  Interpreted  as  being  the  absolute 
value • 
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ABSTRACT 

The  problem  of  deducing  the  values  of  the  exchange  integrals  in  yttrium-  and 
gadolinium-  iron  garnet  from  measurements  of  the  magnetization  and  the  magnetic  con¬ 
tribution  to  the  specific  heat  at  low  temperatures  is  considered.  For  these  garnets  the  spin 
wave  normal  modes  can  be  found  by  solving  the  semi -classical  equations  of  motion  which 
give  rise  to  a  set  of  n  simultaneous  linear  equations,  where  n  is  the  number  of  magnet¬ 
ically  inequivalent  ions  in  the  unit  cell.  Expressions  for  the  thermodynamic  functions  at 
low  temperatures  in  terms  of  the  frequencies  of  the  normal  modes  are  given  assuming  the 
validity  of  the  spin  wave  approximation.  It  is  argued  that  the  temperature  variation  of 
the  frequency  of  these  normal  modes  on  the  macroscopic  properties  con  be  completely 
accounted  for  without  considering  the  zero  point  energy  explicitly. 

Due  to  the  size  of  the  unit  cell,  the  equations  for  the  frequencies  of  the  normal 
modes  can  only  be  solved  numerically  for  general  values  of  l<  .  Such  solutions  ore  ob¬ 
tained  for  k  lying  along  a  [ill]  direction  for  various  values  of  the  exchange  integrals, 
and  fhe  thermodynamic  functions  corresponding  to  these  choices  of  parameters  are 

calculated.  In  the  cose  of  yttrium  iron  garnet,  the  value  of  D,  the  coefficient  of 

2  2 

a  k  in  the  acoustic  dispersion  law,  is  reliably  known  and  fixes  one  linear  combination 

of  J  ,  J  1  /  and  J  , ,  .  By  comporing  our  calculations  with  the  magnetization  data  of 
aa  ad  dd 

Solt,  it  was  established  that  but  since  the  magnetization  was  not  very 

sensitive  to  variations  of  the  ratio  value  could  not  be  estimated  precisely. 

,  -1  -1 
Taking  gives  =  6.35  cm  and  *  31 .8  cm 

o 

For  GdiG  the  specific  heat  data  below  20  K  is  not  very  much  influenced  by  the 
exact  values  of  the  iron-iron  exchange  integrals  which  were  taken  to  be  those  quoted 
above  for  yttrium  iron  garnet.  Again  one  combination  of  and  is  known  from 
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th«  calorimetric  determination  of  the  single  ion  splitting.  By  comparing  the  specific  heat 

data  below  5°  K  with  calculotions  for  various  values  of  J  /J,  it  was  possible  to  determine 

ac  dc  ^ 

J  and  Jj  separately:  J.  *  7. 00  cm  end  J  =  i  75  cm  .  These  values  are  about 

ac  dc  dc  ac 

25%  larger  than  what  one  would  expect  using  the  Weiss  molecular  field  approximation. 
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I.  INTRODUCTION 


Recently  many  investigations  of  the  behavior  of  the  series  of  iron  garnet  compounds 
i 

hove  been  mode.  Pouthenet  first  interpreted  their  magnetic  properties  below  the  Neel 

point  at  550^K  using  the  Weiss  molecular  field  (WMF)  approximation.  More  recently  it  has 
2-5 

been  demonstrated  that  at  the  lowest  temperatures  one  must  use  the  spin  wave  approxi¬ 
mation  in  many  cases  to  interpret  the  thermal  and  magnetic  properties  of  the  iron  garnets. 

3,4 

The  earliest  calculations  gave  a  value  of  the  excitation  energy  of  the  various  normal 

modes  (spin  waves)  for  *  0,  where  I*  a  vector  of  the  first  Brillouin  zone.  While  this 

is  often  sufficient  to  describe  the  resonance  behavior  of  the  iron  garnets,  it  would  seem 

desirable  to  moke  more  accurate  calculations  of  the  macroscopic  properties  which  take  into 

5 

account  the  dependence  of  the  excitation  energies.  Tinkham  has  made  such  a  calcu¬ 
lation  for  a  simplified  model  of  the  interactions  in  ytterbium  iron  garnet  (YbIG).  He  has 
shown  that  there  are  spin  wave  modes  whose  energy  Is  not  very  dependent  on  1^  and  is 

roughly  equal  to  the  energy  of  the  single  ion  splitting  of  the  rare  earth  ion  in  the  WMF 

,,  6 

approximation,  which  can  be  determined  calorimetricoiiy .  Thus  the  physical  picture  of 
the  magnetic  behavior  of  the  iron  garnets  can  be  considered  to  be  well  understood.  The 
aim  of  the  present  more  detailed  calculation  is  twofold.  First,  it  was  felt  that  by  making 
accurate  calculations  of  the  spin  wave  spectrum  for  0  a  more  detailed  comparison 
between  theory  and  experiment  could  be  made.  Second,  it  was  of  interest  to  see  the 
correspondence  in  detail  between  such  a  complete  calculation  and  other  calculations,  such 
as  that  using  the  WMF  approximation,  which  do  not  take  full  cognizance  of  the  structure 
of  the  unit  cell . 


For  yttrium  iron  garnet  (YIG)  and  gadolinium  iron  garnet  (GdIG)  the  difficulty 
in  such  a  calculation  arises  from  the  complexity  of  the  unit  cell .  As  is  well  known,  the 
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frequencies  of  fhe  normal  modes  of  on  isotropic  spin  system  can  be  found  by  solving  the 
classical  equations  of  motion,  which  give  rise  to  a  secular  equation  whose  degree  is  in 
general  equal  to  the  number  of  magnetic  ions  in  the  unit  cell .  A  considerable  simplifi¬ 
cation  in  the  numerical  calculations  is  obtained  if  one  assumes  that  the  energy  surfaces 
in  k  space  are  spherical .  Under  this  assumption  it  is  only  neceuary  to  solve  the  secular 
equation  for  k  lying  along  a  [111]  direction,  in  which  case  the  secular  equation  can 
be  factored.  It  seems  unlikely  that  this  simplification  could  introduce  significant  errors 
into  the  calculation  of  the  macroscopic  properties,  because  the  thermodynamic  functions 
do  not  depend  sensitively  on  the  exact  details  of  the  frequency  spectrum. 

II .  CALCULATION  OF  THE  SPIN  WAVE  SPECTRUM 
A.  Lattice  Structure 

Before  calculating  the  excitation  spectrum  of  spin  waves,  it  is  necessary  to  dis¬ 
cuss  briefly  the  lattice  structure  of  the  iron  garnets.  The  crystal  structure  of  the  iron 

10  7 

garnets  is  cubic,  the  space  group  being  -Ia3d.  The  most  important  symmetry  pro¬ 
perty  of  the  crystal  is  that  the  [ill]  direction  is  a  threefold  axis,  in  Table  I  we  give 
the  positions  of  the  magnetic  ions  in  the  unit  cell .  Each  unit  cell  contains  four  formula 
units  of  iron  garnet,  5Fe202-3M202 ,  where  M  i,s  yttrium  or  any  of  the  rare  earths  from 
samarium  to  lutetium.  Studies  of  the  magnetic  behavior  of  the  garnets  show  thpt  the 

iron  ions  on  the  a  and  d  sublattices  are  strongly  coupled  together  anti ferromagneti cal ly 
/  o 

with  a  resulting  Neel  temperature  of  about  550  K.  When  rare  earth  ions  occupy  the  c 
sites,  their  spin  moments  are  coupled  antiferromagnetically  to  the  resultant  moment  of 
the  a  and  d  sublattices.  This  coupling,  which  is  much  weaker  and  hence  does  not  affect 
the  Neel  temperature,  produces  anomalies  in  the  specific  heat  and  magnetization  below 


-4- 


30°K.  The  magnetic  ions  are  also  si^jected  to  a  crystalline  eiectric  field,  which  often 

8,9 

produces  an  easy  axis  of  magnetization  along  a  [111]  direction.  For  ions  which  are 

3+  3+ 

in  an  S  state,  e.g.,  for  Fe  or  Gd  ,  the  strength  of  this  field  corresponds  to  splittings 
-1  8 

of  the  order  of  .01  cm  or  less  and  hence  can  be  neglected  for  the  present  calculations. 


B.  The  Equations  of  Motion  of  an  Isotropic  Spin  System 

We  now  treat  the  cose  of  an  isotropic  system  of  interacting  spins  in  an  external 
magnetic  field,  H  ,  for  which  the  Hamiltonian  is. 


-s,, -rs,  ^H-s, 

r  r  r  ^ 


(1) 


where  S  and  g  are  the  spin  and  g  value  of  the  ion  at  r  ,  /9  is  the  Bohr  magneton, 

«ke 

and  J  ,  is  an  interaction  coefficient.  The  ions  of  the  three  magnetic  sublattices  are 
r  r 

assumed  to  interact  with  one  another  only  if  they  are  the  closest  pair  in  the  sublattices  in 
question.  The  dipolar  interactions  correspond  to  energies  which  are  very  much  smaller 
than  iron-rare  earth  exchange  energy  and  hence  are  taken  into  occount  only  insofar  as 
their  effect  is  equivalent  to  that  of  a  demagnetizing  field  which  can  be  included  in  H . 
Semi -classically  one  uses  the  torque  equation  to  determine  the  frequencies  of  the  normal 
modes: 


dS 

••r 


=  y  (s  X  H  ) 

- —  r  ••r  ••r 

t  m,  mm  ^ 


(2a) 


dS 


(2b) 


where  y  is  the  gyromagnetic  ratio  of  the  ion  at  r  .  Quantum  mechanically,  the  equations 
r 
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of  motion  ore 


<n|CK  ,  ]|m>  =  Kto  <n|  Sj.  |m> 


(3a) 


and 


<n|CK,  ]|m>  *  -Kw  <n|  |m>  , 

1  X  ,  y 

wh«r«  ±  iS^  .  One  finds,  using  Eq.  1  and  the  commutation  relations. 


S  X  S  =  i  S  ,  that 

Ota 


[K,S  '^]  =  g  j9S  -  2r  J  ,(S^^S  -  S^S  . 

r  r  r'  JL  ’’  «L 


and 


(3b) 


(4a) 


z 


C3C,S  ]  »  -g  +  2r  J  ,(S  S  ,  -S  S  , 
r  '’r  r  r  r  r  r  r  r 


)  , 


(4b) 


if  the  magnetic  field  is  oriented  in  the  plus  z  direction.  In  the  random  phase  approxi¬ 
mation  (RPA)  one  replaces  the  operator  by  its  thermal  value,  (S^  ).  This 
corresponds  exactly  to  linearizing  the  semiclossical  equations  of  motion. 

In  either  cose  if  one  writes  S,  =(2Nyc)  ^  T  S  e  "*  ^  , 

k  T  R  R  +T 

^  0^ 

where  t  denotes  the  position  of  an  ion  in  the  unit  cell,  ^  is  a  translation  vector  of 

the  lattice  (which  we  take  to  be  body-centered  cubic),  |<  is  a  vector  in  the  first 

Brillouin  zone,  and  N  is  the  number  of  unit  cells,  then  one  finds^  that  for  each 
uc 

value  of  l<  the  allowed  values  of  Koj  are  the  eigenvalues  of  the  matrix  A ,  whose 
components  are 

- 


A.j  =  Cg,/5H-2S  y...  (0)]  6,. +2 J,.  (^)  S.  yii(k) 


(5) 


ik  •  (r  -  r.,) 

where  the  indices  i  and  j  label  the  ions  in  the  unit  cell,  and  y.,  =  e  **  '  **1  , 

the  sum  being  taken  over  nearest  neighboring  ions  in  the  sublatticf .  The  dimensionality 
of  this  matrix  is  equal  to  the  number  of  magnetic  sublattices,  i  .e.,  20  for  YIG  and  32  for 
GdiG.  These  matrices  are  shown  in  Fig.  1  and  2,  where  the  submatrices  are  given  in  Fig. 
3-7.  The  numbering  of  the  rows  and  columns  of  A.,  correspond  to  the  numbering  of  the 


•I 


magnetic  ions  as  given  in  Table  I . 


Let  us  consider  some  generol  properties  of  this  matrix.  A,,  .  It  can  be  transformed 


'I 


into  the  form. 


B  C 
-Z*  D 


(6) 


where  D  and  ]3  are  Hermition  submatrices.  To  diagonalize  A ,  we  look  for  a  matrix  0 
such  that 


O  AO  *  A,  , 

m,  mm  •.() 


(7) 


where  A^  is  diagonal,  and  where,  due  to  the  form  of  A ,  O  and  O  satisfy  the  con¬ 
ditions 


0  0  =  1  ,  O  =  ( 

ehe  eta  mm 


-11  2i2 
^21  ^22 


6*  -  3* 

t  xn  X21 

)  o  =  (  -  ~  )  . 

'  -O*  o*  ' 

-12  22 


(8) 


where  the  O..  are  suitable  submatrices.  Since  A  is  not  Hermition,  the  roots  need  not 

-II 

3 

be  real,  but,  as  indicated  by  Dougioss,  this  would  imply  an  unstable  ground  configuration. 

In  terms  of  the  occupation  numbers,  n  ,  ,  the  energy  levels  of  the  spin  system  are  ap- 

Pi 

proximately 


I 

I 


% 
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00 

where  p  labels  the  branch  of  the  spectrum  and  E  is  the  ground  state  energy.  In  the 

Appendix  we  show  that  the  temperature  dependence  of  the  frequencies  of  the  normal  modes 

12 

does  not  affect  the  zero  point  energy. 

Let  us  now  calculate  u  ,  the  change  in  the  magnetic  moment  when  n  ■  is 
PK  Pa 

increased  by  one,  remembering  that  the  magnetic  moment  of  a  state  is  given  by  the  derivative 

13 

with  respect  to  the  magnetic  field  of  the  energy  of  that  state.  Then 

d 


"ok  '  dH“",k^'>'*'"okl-'okl>'“„k 

••  ^  ^ 


(lOo) 


■^IX“okl*%-^*<"„k  ' 


(10b) 


where  a  ,  the  polarization  of  the  mode,  is  just  the  sign  of  Kw.i,  •  Lking  pertur- 
p  PJi 

bation  theory  one  finds 


f  dA(k)  t 

(O  O)  =  a  j8(0  GO)  , 

““dH  Po  o~“''^PP 


(11) 


where  G  is  diagonal  and  G  is  the  g  value  of  the  m*^^  ion  in  the  unit  cell .  When 

mm 

all  the  g  values  are  the  same,  as  for  YIG  and  GdIG,  one  has  the  simple  result  that 


M  =  P  j8  .  (12) 

®h,  ® 

We  now  consider  some  properties  of  the  solution  to  this  eigenvalue  problem  taking 

account  of  the  symmetry  of  the  gcrnet  lattice.  For  1^  =0 ,  the  symmetry  is  sufficient 

to  permit  a  factorization  of  the  secular  determinant  and  thence  a  complete  analytic 

determination  when  k  lies  along  a  [1 1 1]  direction,  which  we  assume  to  be  the  case  in 

the  following  discussion.  We  denote  by  R  the  operator  which  rotates  the  crystal  about 

2ff 

the  [ill]  direction  by  “J"  •  We  note  that  since  commutes  with  A,  the  eigenvectors 
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of  A  are  also  eigenvectors  of  .  Since  =  1  ,  the  eigenvalues  of  R  are  just  the  cube 

2  2vi 

roofs  of  unity,  1,  X ,  and  X  ,  where  X  =  e  — g —  .  The  secular  equation  may  then  be 
written  in  the  form. 


Det|A(j<,  1)-Kw^  ,l|A(k,X)-Kw^  ^^llA(k,X  ^2  \  -  O 


(13) 


where  the  eigenvectors  of  A(l(  ,  r)  are  associated  with  the  eigenvalue  r  of  R  .  These 
matrices  are  given  in  Fig.  8  and  9,  oikI  in  Table  II  we  give  the  linear  combinations  of 
rows  and  columns  of  A  which  correspond  to  the  rows  and  columns  of  A(|^,  r).  Using 
Fig.  9,  one  verifies  that 

p'’A(k  ,X)P  =  A(-k,X)  =  A*(k,X^)  ,  (14) 


where  P  is  a  unitary  transformation  which  interchanges  rows  and  columns  as  follows: 

3  **  6,  4  s*  5,  7  **  10,  8  **  9.  Since  the  frequencies  are  assumed  to  be  real, 

2 

A(I^,X)  and  A(Jc,X  )  have  the  same  eigenvalues. 

Although  the  secular  equation  factors,  analytic  solutions  are  obtoinoble  only 

3,4 

at  the  center  or  extreme  corner  of  the  Brillouin  zone.  However,  since  the  mag¬ 
netic  properties  at  the  lowest  temperatures  depend  only  on  the  behavior  of  the  low 
energy  part  of  the  spectrum  it  may  be  useful  to  exapnd  the  energy  of  the  acoustic  branch 
in  powers  of  k.  Due  to  the  cubic  symmetry  this  expansion  takes  the  form 


22  A4  4  22  22  22 

Kco=gj8H  +  Da  k  +Ea  k  +Fa  (k  k  +k  k  +k  k  )  + 

X  y  X  z  y  2 


(15) 


By  perturbation  methods  one  finds 
D 


200J  -  125J  . +75J..  -  20J.  S  +  50J  S 

aa  ad  dd  dc  c  ac  c 


16(  -5  +  6S^) 


C 


(16) 
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For  YIG  (S  =0)  one  also  finds 
c 


12288 


C(32J«,  -  19Jod  +  - 


(48J 


ao 


96  J 


F  -  12288 


8J  -  6J  . 
oo  ad 


^Jgd  ^ 
Jod 

] 


(17) 

(18) 


For  GdIG  the  expressions  for  E  and  F  are  somewhat  unwieldy  and  hence  it  may  be 
more  revealing  to  give  approximate  expressions: 


-6S-D  5 

E» - ^ -  +  - C-160Jaa  +  9SJad  -  60Jdd  + 


A(-5  +  6S)  12288(-5  +  6S) 

c  c 


5(48J  -  25J^  +  12Jdd)" 

+ - ^-j - ] 

•’ad 


240J 


F  = 


12288  (-5+6$^) 


C-640J  -  80J  + 


od 


aa  aa  3J  ,  -  4J 
od  aa 


■]  , 


(19) 

(20) 


where  t  =  ^•’ac  '*  energy  level  splitting  of  the  rare  earth  ion 

in  the  WMF  approximation.  For  GdIG  the  term  in  Eq.  15  involving  F  is  much  less 
important  than  that  involving  E  ,  which  can  usually  be  adequately  approximated  by 
the  first  term  of  Eq.  19. 

Several  observations  about  these  expansions  should  be  made.  One  sees  that 
the  first  anisotropic  term  is  that  involving  F,  which  vanishes  for  =  0. 

In  any  event,  one  sees  that  rotio  f/E  is  small  for  reasonable  values  of  the  other 
exchange  integrals,  in  agreement  with  our  assumption  that  the  energy  surfaces  in  k 
space  are  nearly  spherical .  One  also  notes  that  E  has  the  opposite  sign  from  D 
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since  all  the  exchange  integrals  ore  taken  to  be  negative,  with  the  largest. 

It  is  interesting  to  note  the  role  played  by  the  spin  of  the  rare  earth  ions,  •  it  enters 

all  the  expressions  in  the  denominator  via  a  factor  (-5  4  6S  ),  thus  profoundly  affect- 

c 

ing  D  and  hence  the  low  temperature  properties.  The  condition  for  ant! ferromagnetism 
is  just  -5  4  6Sg  ==  0,  in  which  case  our  expreuions  are  no  longer  valid. 

In  order  to  discuss  the  magnetic  properties  at  higher  temperature,  it  is  neceuary 
to  solve  for  the  normal  modes  numerically,  which  necessitated  the  use  of  an  electronic 
computer,  considering  the  size  of  the  matrices  involved.  We  used  a  method  of  succes¬ 
sive  approximations:  we  eliminated  the  largest  off-diagonal  element,  ,  by  per¬ 
forming  successive  two  dimensional  rotations,  i.e.,  by  performing  transformations  with 
a  matrix  O  obeying  Eq.  iO  and  of  the  form. 


0..  «  fij. 

Il  '1 

i  /  n,m 

1  ^  n,m 

(2ia) 

O  .  =  0 

1^  1 

i  n,m 

0  .  =  0 
mi 

1  /  n,m  . 

(21b) 

The  numerical  errors  incurred  in  these  calculations  were  negligibie,  as  was  verified 

by  comparing  the  trace  of  the  original  matrix  with  the  sum  of  the  eigenvalues.  In 

addition,  for  |<  =  0 ,  we  found  agreement  between  our  numerical  results  and  the 

4 

analytic  expressions  of  Dreyfus. 

The  results  of  calculations  for  various  values  of  the  ratios  o''®* 

Jdd/Jod  of  YIG  are  given  in  Fig.  10,  11,  and  12.  In  Section  IV  we 

discuss  the  determination  of  for  given  values  of  these  ratios.  The  most  striking 
feature  of  the  spectrum  is  that  the  frequency  of  most  of  the  modes  does  not  depend 
strongly  on  k  .  We  can  compare  the  frequencies  with  those  one  would  expect  using 


Ka//y 


« 


(cm"') 


-13- 


the  WMF  model .  Here  one  would  hove  eight  Iron  ions  on  a  sites  in  on  effective  field, 

H  ,  where 

Q 

9^^  =  -30Jad^^-'oo  (22a) 

and  12  iron  ions  on  d  sites  in  on  effective  field,  ,  where 

From  the  grophs  one  sees  that  on  the  whole  the  frequencies  lie  somewhat  below  those  of 
the  WMF  model  and  hence  the  values  of  the  exchdnge  integrals  os  deduced  from  com¬ 
parison  of  experimental  data  with  the  WMF  theory  are  expected  to  be  less  than  those  we 
find.  This  discrepancy  is  the  most  serious  at  the  lowest  temperatures,  just  where  precise 
magnetization  data  is  availdble. 

In  Fig.  13  and  14  we  show  the  low  frequency  part  of  the  spectrum  for  GdiG  for 

k  lying  along  a  [ill]  direction.  These  low  optical  modes  correspond  in  the  WMF 

3+ 

picture  to  the  reversal  of  o  Gd  spin  in  the  effective  magnetic  field  of  the  iron  ions. 

14 

The  high  frequency  part  of  the  spectrum  is  not  very  different  from  that  of  YIG.  We 

15 

have  selected  values  of  and  such  that  the  average  frequency  of  the  low 

6 

optical  modes  is  that  determined  calorimetricaily .  Fig.  13  and  14  therefore  show  the 
effect  on  the  spectrum  of  varying  the  ratio  keeping  the  average  of  the  low 

optical  frequencies  constant.  For  a  given  vaiue  of  this  average  frequency  varying  the 
iron-iron  exchange  integrals  has  relatively  little  effect  on  the  low  frequency  part  of 
the  spectrum.  In  contrast,  as  the  ratio  increased,  the  low  optical  modes 

are  split  opart  and  the  average  frequency  becomes  significantly  less  than  the  WMF  split¬ 
ting  which  is  given  by 

gjSH  =-10J.  ^20J 
®  c  dc  oc 


(23) 


-14- 


These  results  agree  with  the  analytic  expressions  found  by  Dreyfus  for  k  =  0,  which  show 

2 

that  the  spread  in  the  low  optical  modes  is  roughly  proportional  to  where  is 

an  iron -rare  earth  exchange  integral.' 

III.  CALCULATION  OF  THE  MACROSCOPIC  PROPERTIES 


Once  the  spin  wave  spectrum  has  been  found,  the  magnetization  and  the  specific 
16 

heat  are  easily  calculated.  For  instance,  for  the  magnetization  one  has 


M(t)  =  M(0)  -  — Y’  ^ 
Sir'*  n 


; 

B.Z. 


1 

1  -exp(-x„(k)) 


(24) 


where  the  sum  is  taken  over  the  normal  modes,  n ,  the  integration  over  the  Brillouin  zone, 

Ka)-(k) 

fi  (k)  is  given  by  Eq.  11  and  x  (k)  = - ^  where  k.  is  the  Boltzmann  constant. 

“b' 

In  order  to  calculate  M(T)  from  the  spectrum  os  determined  for  j<  lying  along  a  [111] 
direction,  we  assumed  spherical  energy  surfaces  and  integrate  over  a  sphere  of  radius 
k 


4irk  ^  2ir  3 

such  that  "’OX  =  2( - )  since  there  are  two  ions  per  unit  cell .  Therefore 

max  2  o 


max 


1 


M(T)=M(0)--^r  f  •  PS 


For  YIG  and  GdIG  where  all  the  ions  have  the  some  g  value,  one  can  calculate  the 
magnetization  at  low  temperature  when  the  expansion  of  Eq.  15  is  valid: 

M(0)  kftT  3/2  3  .  kftT  5E+F 

AM  =  M{0)-M(T)=  CC (3/2) -■7C(5/2)-l- (-ir“)]  ,  (26) 

^TOT  ^  D 


where  is  the  totai  spin  per  unit  cell,  i«.,  *  20  for  YIG  and  * 

for  GdIG.  In  an  analogous  way  one  finds  the  specific  heat  to  be 


32 
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C  =  k 


B 


V 


'^max  2  , 

E  J*  [x  (k)3  - 2  4ffk^dk  , 

"  0  "  [1  -exp(-x^(k))r 


(27) 


v/here  V  is  the  volume.  Again  in  zero  applied  field  one  can  calculate  the  specific  heat 
when  the  expansion  of  Eq.  15  is  valid; 


^  ^  CC  (5/2)  -  7/4:(7/2)  ^ 

Dir 


(28) 


Expansions  for  the  case  when  there  is  an  applied  field  have  been  given  by  Robinson^^ 
assuming  a  quadratic  dispersion  relation.  From  his  graphs  one  can  see  the  relatively 
large  effect  of  a  magnetic  field  on  the  magnetization. 

It  is  interesting  to  consider  the  range  of  validity  of  the  expressions  (26)  and  (28). 
Taking  ■^dc/”*ad  *  rough  accord  with  experi¬ 

mental  data,  one  finds  using  Eq.  16,  17,  and  19  that  E/D  -  .31  for  GdIG  and  .12 

3/2 

for  YIG.  Using  Eq.  27  we  find  a  10%  deviation  from  pure  T  behavior  for  the  specific 

heat  when  kg  T  w  D/8  for  YIG  and  kgT  «  .04 D  for  GdIG.  Since  DAg  45®  K  for 

o  3/2 

YIG  and  15  K  for  GdIG,  one  sees  that  whereas  YIG  exhibits  T  behavior  below  say 

3/2  o 

6°K,  GdIG  may  never  really  obey  a  T  ^  law  since  below  1  K  anisotropy  ond  dipolar 

18 

perturbations  will  influence  the  thermal  properties. 

We  now  estimate  the  effect  of  the  temperature  dependence  on  our  calculations. 

In  the  Appendix  we  show  that  at  low  temperatures  one  should  use  the  average  frequencies, 
Ka),  to  calculate  the  thermodynamic  functions  where 


^  Kft)(0)  •^Kctf(T) 

Kw  (T)  *  2 


(29) 


To  calculate  )(w(T)  one  simply  replaces  S|  by  its  thermal  value;  however,  for  the  iron 
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ions  below  soy  SO^’k  the  thermal  value  does  not  differ  appreciably  from  the  value  at 

T  =  O^’k.  In  Table  Hi  we  compare  the  frequencies  of  the  normal  modes,  ji(u,  calculated 

for  T  =  0®  K  with  those  calculated  fe  T  =  20®  K.  The  small  temperature  variation  of 

those  frequencies  is  not  experimentally  detectible  at  present.  It  is  interesting  to  note 

the  increase  in  the  frequencies  of  the  low  optical  modes  as  the  temperature  is  raised  and 

the  RPA  and  the  WMF  become  equivalent.  This  effect  results  from  the  change  in  the 

effective  field  acting  on  the  iron  ions  due  to  the  slight  temperature  dependence  of  the 

rare  earth  magnetization  below  20°  K,  and  is  not  to  be  confused  with  the  effect  of  rare 

earth-rare  earth  interactions,' or  with  the  effect  of  the  temperature  dependence  of  the 
iron  sublattice  magnetization  both  of  which  we  have  neglected. 

IV.  THE  DETERMINATION  OF  THE  EXCHANGE  INTEGRALS 
FROM  EXPERIMENTAL  DATA 

1.  YIG 


Yttrium  iron  garnet  has  been  the  object  of  several  experimental  studies  from 

which  information  about  the  exchange  integrals  could  be  obtained.  Using  o  molecular 

1 

field  analysis,  Pauthenet  was  able  to  determine  values  of  the  exchange  integrals  which 

19 

were  later  apparently  confirmed  by  high  temperature  susceptibility  measurements. 

However,  the  values  of  the  exchange  integrals  so  determined  give  a  value  of  15  cm 

for  D  (see  Eq.  15)  which  is  in  disagreement  with  several  subsequent  determination  of 

this  quantity.  For  instance,  by  comparing  the  experimentally  determined  low  tempera- 

6,  20-22 

ture  specific  heat  with  the  results  of  the  spin  wave  calculation,  one  was  able 

to  deduce  that  D  «  27  cm  .  This  value  of  D  has  been  corroborated  by  observations 

23  24 

of  the  microwave  instability  in  YIG.  '  Since  the  determination  of  D  using  a  spin 


wave  theory  i.s  the  more  unambiguous  both  from  a  theoretical  and  on  experimental 


-17- 


standpoint,  we  assume  it  to  be  the  more  reliable.  However,  fixing  the  value  of  D  to  be 
27  cm  ^  does  not  serve  to  determine  the  exchange  integrals  uniquely  since  from  Eq  .  15 
we  have: 


D  =  -~(8J  -5J  .  +  3J,.)  . 

15  aa  ad  dd 


(30) 


The  problem  we  consider  is  to  determine  values  for  ,  and  consistent  with 

the  known  value  of  D  and  which  best  reproduce  the  detailed  behavior  of  YIG. 

25 

Recently  Wojtowicz  has  suggested  that  by  taking  J  =  J  ,  ,  =  0,  one  might 

QQ  do 

be  able  to  fit  the  susceptibility  data  cbove  the  Neel  point  using  a  linked  cluster  ex¬ 
pansion  rather  than  a  molecular  field  approximation.  He  was  able  to  interpret  the  ex¬ 
perimental  data  in  this  way,  but  using  a  value  of  which  corresponds  to  D  =  38  cm 
It  is  quite  possible,  however,  that  one  could  also  fit  the  data  using  non-zero  values  of 
and  ,  although  the  analysis  would  be  formidable.  In  any  event,  exact  agree¬ 
ment  between  determinations  of  D  over  widely  separated  temperature  intervals  is  not  to 
be  expected,  due  to  the  sensitive  dependence  of  the  exchange  integrals  on  the  lattice 
constant.  Also,  the  variation  of  D  with  temperature,  as  determined  by  microwave  in- 
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stability  measurements,  does  not  seem  to  be  consistent  with  =  Jjj  =  0 ,  at  least 

.  .  24 

according  to  the  RPA. 

26 

Recently  I .  Solt  has  made  accurate  measurements  of  the  magnetization  of  YIG 

in  a  magnetic  field  of  4000G  for  temperatures  between  5  and  50°  K  using  the  properties 

of  magnetostatic  modes.  Since  in  this  temperature  range  the  magnetization  varies  by 

only  .4%,  one  should  check  that  volume  changes  do  not  influence  the  magnetization 

27 

appreciably.  In  fact,  using  the  GrUneisen  relation  and  taking  the  compressibility  to 

28 

be  roughly  that  at  room  temperature  os  measured  by  Kaminow,  one  finds  the  change  in 
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the  magnetization  due  to  the  explicit  volume  dependence  to  be 


iiM 

M 


ty  _ 


5C 


(3)(7x  iC"'^ 


=  ^  ]■  C 


Lattice 


dT  « 


3C0  cm 


—  (3^;  joules)^  2^  X  10  , 

(31) 


where  y  is  the  GrOn^isen  consiant  end  x  the  compressibility,  so  that  the  volume  change 

is  responsible  for  only  i%  of  iha  observed  decrease  in  the  magnetization  at  5t°K.  (The 

lattice  specific  heat  was  estimated  from  the  data  of  Harris  and  Meyer^) .  It  is  also  known 

that  the  temperature  dependence  of  the  exchange  integrals  due  to  the  change  in  the  lat- 

29 

tice  constant  in  this  temperature  range  is  negligible. 

However,  the  application  of  an  external  field  of  even  4lCt  G  affects  the  mag- 
netiswtion  signifieently  in  this  temperatures,  os  we  have  mentioned  previously.  The 
magnetization  was  measured  by  measuring  the  fields  for  resonant  excitation  of  the  (210) 

and  (220)  magnetostatic  modes  at  a  frequency  of  9092  me.  The  resonance  conditions 

26 

are 


8  ir 

-  »/y  -  . 


(32a) 

(32b) 


where  M(H)  is  the  magnetization  M  in  a  field  H  and  H  is  the  field  for  resonance 

Imn 

of  the  (Imn)  magnetostatic  mode.  The  magnetization  found  by  subtracting  the  second 
equation  from  the  first  corresponds  to  an  overage  field  H  for  which 

*^210  "^220  "  1T^^^210^'^15^^”220^  “ 
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Expandlng  M(H2jq)  and  about  H,  one  finds 


H  =  220  210  =  ^/y  +  M 

3 


(34) 


30 

Anderson  and  Suhl  have  found  a  dispersion  relotion  for  the  acoustic  mode  taking 
account  of  dipolar  interaction,  which  for  a  sphere  can  be  written, 

(K<*>)  =  CDa\^ +g/3(H  -  M)]  [Da^  +gj9(H  -  ~:^M +4irM  sin^  0/  )]  , 


(35) 


where  is  the  angle  between  j<  and  the  magnetic  field,  Since 
4irM/3H  .2  one  finds 

Kw  «  Da^k^  +  g/J(H  -  M  +  2ir  M  sin^  0-|^  )  . 


(36) 


For  a  given  value  of  |  k  |  this  dispersion  relation  introduces  a  breadth  of  2irg/3M 

into  the  spectrum,  but  leaves  the  average  frequency  unchanged.  Accordingly,  the 

demagnetizing  field  can  be  neglected  in  the  first  approximation.  In  addition  we  have 

assumed  that  the  effects  of  anisotropy  can  be  taken  account  of  by  the  usual  anisotropy 
31 

field,  -4K/3M,  where  K  is  the  first  order  onisotropy  constant.  Taking 

8 

4irM/5  =  487Q  and  K/M  =  130G  we  find  the  total  field  to  be: 


4K  -1 

—  =  .302  cm 


(37) 


In  order  to  determine  AM(T)  =  M(0)  -  M(T)  it  is  necessary  to  extrapolate 
M(T)  to  T  =  0 ,  which,  however,  can  not  be  done  without  introducing  a  significant 
error.  In  the  low  temperature  limit  M(T)  depends  only  on  ^  and  D,  which  are  known, 
so  that  M(0)  could  be  determined  by  comparing  the  values  of  AM  for  temperatures  below 
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10  K  with  the  calculated  values.  The  uncertainty  in  this  procedure  did  not  affect  our 
determinations  of  the  exchange  integrals.  We  calculated  the  magnetization  for  various 
values  of  the  ratios  J ond  the  value  of  being  determined  to  give 

the  known  value  of  D.  By  comparing  the  family  of  curves  so  obtained  with  the  ex¬ 
perimental  data  it  was  then  hoped  to  determine  J  and  J,,.  It  happened,  however, 

aa  dd 

that  the  magnetization  was  not  sensitive  to  variations  of  the  rotio  •  O'’ 

other  hand  the  value  of  J  /J  .  did  influence  the  calculated  magnetization  curves, 

aa  od 

os  can  be  seen  in  Fig.  15.  We  were  thus  able  to  deduce  J  /J  ,  ■  .2,  in  contrast 

oa  od 
25 

to  the  approximate  analysis  of  Wo{towicz.  Considering  the  geometry  of  the  garnet 
lottice,  it  would  be  surprising  if  the  exchange  coupling  between  neighboring  d  ions 
were  weaker  than  that  between  neighboring  a  ions,  since  the  latter  are  farther  apart 
than  the  former,  so  that  we  assume  ■•dd  r  -*00  ^ .  Since  there  is  an  optical  modes  whose 
frequency  for  ^  ■  0  Is  20  -  40  wo  were  able  to  estdblish  the  upper  bound 

Jjj  <  «4J^  •  In  Table  IV  we  give  the  corresponding  sets  of  exchange  integrals  together 
with  those  of  GdIG  whose  determination  is  discussed  below.  It  is  hoped  in  the  future  to 
estimate  the  ratio  known^  value  of  the  Ndel  temperature. 

2.  GdIG 


As  we  have  mentioned  previously,  the  details  of  the  low  frequency  part  of  the 
spectrum  for  GdIG  are  not  very  sensitive  to  the  exact  values  of  the  iron-iron  exchange 
integrals  (assuming  D  to  be  knoWh).  We  therefore,  somewhat  arbitrarily,  mode  calcu¬ 
lations  taking  “  *2  and  ~  *2  In  the  first  case  and  Jjy/Jgj  “  *3 

in  the  other.  The  object  of  such  calculations  was  to  determine  and  separately 


by  fitting  the  experimental  results.  Dreyfus  has  previously  determined  values  of  J 


ac 


lO^x  AM/M 
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and  using  the  frequencies  for  .1^  «  0  for  a  particular  value  of  the  ratio  J  • 

Although  this  procedure  does  not  cause  a  great  error  in  the  determination  of  the  exchange 

integrals,  we  thought  it  worthwhile  to  attempt  such  a  determination  taking  account  of 

the  k-dependence  of  the  low  optical  modes  and  also  seeing  to  what  extent  varying  the 

ratio  J  /J  influenced  the  specific  heat.  Above  10  K  the  specific  heat  (per  unit  cell) 
ac  dc 

is  well  approximated  by  24  Einstein  functions  appropriate  to  the  average  frequency  for 

2 

large  values  of  k  ,  since  the  factor  k  dk  in  the  density  of  states  weights  the  large  values 
of  k  the  most.  This  average  frequency  is  known  accurately  from  calorimetric  measure¬ 
ments^  and  for  a  given  value  of  J  /J  ,  fixes  the  values  of  these  exchange  integrals. 

We  then  compared  the  experimental  and  theoretical  values  of  the  specific  heat  below 

5°K  for  various  values  of  J  /j  ,  .  The  results  of  such  a  procedure  are  shown  in  Fig. 

ac  dc 

16.  It  was  found  that  J  /J  ,  *  .25  gave  the  best  fit  to  the  data  for  both  sets  of 
ac  dc 

values  of  the  iron-iron  exchange  integrals  considered.  The  corresponding  values  of 

J  and  J  which  depend  slightly  on  the  choice  of  the  values  of  the  iron-iron  exchange 
ac  dc 

integrals  are  given  in  Table  IV  and  are  about  25%  larger  than  one  would  expect  using 

the  WMF  approximation.  The  reason  for  this  can  be  seen  from  Fig.  13  and  14  where  we 

indicate  both  the  frequency  corresponding  to  the  WMF  acting  on  a  rare  earth  ion  as  given 

by  Eq.  23  and  the  average  frequency  of  all  the  low  optical  modes.  Equating  the  frequency 

of  the  single  ion  splitting  as  determined  calorimetrically  with  the  WMF  frequency  clearly 

leads  to  smaller  values  of  the  exchange  integrals  than  we  find.  One  can  also  compare 

-1 

the  values  of  J  and  J  ,  we  fmd  with  those  found  by  Dreyfus:  J  =0.49  cm  and 
ac  dc  '  '  ac 

J  =  4. 17  cm  .  This  surprising  discrepancy  between  his  values  and  ours  is  mostly  due 
dc 

to  a  numerical  error  in  his  calculation  of  the  lowest  optical  mode  which  unfortunately  ^ 
affects  the  determinotion  of  the  exchange  integrals  rather  critically. 


e 


SPECIFIC  HEAT  (JOULES/ MOLE  ‘K) 
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It  should  be  noted  that  we  were  unable  to  obtain  a  good  fit  to  the  specific  heat 

data  for  temperatures  below  3°K.  The  reason  for  this  discrepancy  is  not  clear  at  present 

and  will  be  investigated  experimentally  in  the  near  future.  We  did  not  attempt  to  obtain 

6 

a  better  fit  by  altering  the  value  of  D  as  we  did  previously  which  accounts  for  most  of 

the  difference  between  the  calculations  presented  here  and  those  given  previously.  Above 

10°K  the  calculated  and  experimental  values  of  the  specific  heat  are  in  close  agreement. 

As  we  show  in  the  Appendix,  one  does  not  expect  that  the  small  temperature  dependence 

of  the  frequencies  of  the  normal  modes  will  influence  the  specific  heat  significantly. 

In  Fig.  17  we  show  the  magnetizotion  of  GdiG  below  20  K  as  calculated  for  one 

of  the  sets  of  exchange  integrals  given  in  Tdble  IV.  The  curve  calculated  for  the  other 

set  of  exchange  integrals  is  indistinguishable  from  that  shown.  The  magnetization  of  GdIG 

1  32 

has  been  measured  by  Pouthenet  and  by  Wolf  and  Bozorth,  but  due  to  experimental 

difficulties  these  measurements  are  not  sufficiently  refined  to  be  suiteble  for  comparison 

with  spin  wove  theory  in  this  temperature  range.  No  information  could  be  obtained  from 

33 

the  absorption  spectrum  of  GdiG  os  measured  by  Sievers  ondjinkhom,  since  the  single 

5 

ion  transitions  ore  forbidden  in  GdiG. 


wv 
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TABLE  i. 

THE  LOCATION  OF  THE  MAGNETIC  IONS  IN 
THE  UNIT  CELL  OF  THE  GARNET  LATTICE® 


a  sites 

d  sites 

c  sites 

1.  0,0,0 

0,1/4, 3/8 

0,1/4, 1/8 

2.  1/2, 1/2,0 

3/8,0, 1/4 

1/8,0, 1/4 

3.  0,1/2, 1/2 

l/4,3/8,0 

l/4J/8,0 

1/2,0, 1/2 

0,3/4, 1/8 

0,3/4,3/8 

5.  1/4, 1/4, 1/4 

1/8,0,3/4 

3/8, 0,3/4 

6.  3/4, 1/4, 3/4 

3/4, 1/8,0 

3/4,3/8,0 

7.  3/4, 3/4, 1/4 

0, 1/4, 7/8 

0,l/4,5/8 

8.  1/4, 3/4, 3/4 

7/8,0, 1/4 

5/8,0, 1/4 

9. 

1/4, 7/8,0 

l/4,5/8,0 

10. 

0,3/4, 5/8 

0,3/4,7/8 

11. 

5/8,0, 3/4 

7/8,0, 3/4 

12. 

3/4,5/8,0 

3/4, 7/8,0 

°  The  coordinates  of  the  magnetic  ions  are  given  os  fractions  of  the  lattice  parameter  which 
is  12.378  %  for  YIG  and  12.465^  for  GdlG.^  The  additional  sites  in  the  unit  cell  are 
obtained  by  adding  (^,  to  the  above  sites.  This  table  is  compiled  from  reference  7. 
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TABLE  II. 


THE  TRANSFORMATION  WHICH  REDUCES  THE  SPIN  WAVE  MATRIX® 


1. 

For  A(l^,  1) 

a. 

For 

2. 

1 

1/^/3  (02  +  ag  4  a^) 

1. 

1/Z3  (02  4  Xog  +  X^a^) 

3. 

®5 

4. 

1//3  (a^  ♦  a^  +  00) 

2. 

1/Z3  (a^  +  Xjp^  4  X^Og) 

5. 

1/^3  (d^  ^2  **3^ 

3. 

l/(</3  (d  ^  +  Xd2  4  X^dg) 

6. 

l/v3  (d^  4  dg  +  d^) 

4. 

1//J  (d^  +  Xd^  4  X^d^) 

7. 

I//3  (d^  4  dg  +  d^) 

5. 

1/73  (d^  4Xdg  +  X%) 

8. 

1/^/3  (djQ  4  d^^  +  d^j) 

6. 

1/73  (r^lQ  4X  d^2) 

9. 

1/^/3  (c^  +  4  Cg) 

7. 

1/73  Cc^  4  XCg  4  X^Cg) 

10. 

V-ZS  (c^  4  Cg  4  c^) 

8. 

1/73  (c^  +XC5  +X^c^) 

11. 

1/Z3  (C^  4  Cg  +  c^) 

9. 

1/73  (C7  +  Xcg  +  X%) 

12. 

l/V3(c,o+c„ 

10. 

‘’The  rows  and  columns  of  the  matrices  A(k,  1)  and  A(k,X)  are  labelled  by  numbere.  The  rows 
and  columns  of  the  original  matrix  are  l*elled  by  “a  letter  and  a  number.  The  letter  tells  the 
sublottice  and  the  numbering  within  the  sublattice  is  according  to  Table  I . 

^’The  linear  combinations  for  ore  found  by  substituting  for  \  everywhere. 


TABLE  III 


TEMPERATURE  DEPENDENCE  OF  THE  FREQUENCIES  OF 


THE  NORMAL  MODES  OF  GdlG® 

r  =  1 

r  =  X 

T-0®K 

T“20®K 

T*0®K 

T=^20®K 

1. 

19.733 

19511 

19J33 

19511 

2. 

9.012 

9.054 

13561 

13534 

3. 

-3.307 

-3.164 

-18741 

-18545 

4. 

-15.083 

-14.978 

-21 .083 

-20.978 

5. 

0.000 

OjOOO 

-12.957 

-12545 

6. 

1.056 

1.057 

-16702 

-16507 

7. 

13.561 

13534 

1.101 

1.101 

8. 

1 1 .687 

1 1 .767 

0759 

0.769 

9. 

-12.957 

-12545 

0720 

0.731 

10. 

-18.741 

-18.646 

1 .056 

1.057 

11. 

0.759 

0.769 

12. 

1,101 

1.101 

®These  frequencies  for  k  =  0,  which  ore  given  in  units  of  |  Jgj  1  ^-are  calculated 
according  to  Eq.  29.  r* is  the  rotation  eigenvalue.  The  modes  with  rotation  eigerfral 
X  are  degenerate  with  those  with  eigenvolue  X  ,  (see  Eq.  13). 


TABLE  IV. 


POSSIBLE  VALUES  OF  THE 
EXCHANGE  INTEGRALS  IN  YIG  AND  GdIG® 


'oa/Jod 

■’dd/'^od 

Jod 

J 

aa 

■•dd 

Jdc 

J 

ac 

.2 

.2 

31.8 

6.35 

6.35 

7.00 

1.75 

.2 

.3 

35.6 

7.12 

10.68 

7.12 

178 

.2 

.4 

40.5 

8.09 

16.18 

not  calculated 

^Values  of  the  exchange  integrals  are  given  in  cm 
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APPENDIX 


Calculation  of  the  Zero  Point  Energy 


34 

To  calculate  the  zero  point  energy  we  use  the  Holstein-Primakoff  expansions 


S,  »  S.  -  a.  a. 

■  III 


S,  =  a/2S  a."  .. 


S.  =  a.  .. . 

I  I 


(AJa) 

(Alb) 

(Ale) 


for  spins  oriented  in  the  plus  z  direction,  and 


z  +  - 

S.  =  -  S.  "t  a.  a. 

I  III 

S  =  ^S.  a  *  .. . 

I  I  I 


(A2a) 

(A2b) 


S.  =  v(2S.  a  ... 
I  i  < 


(A2c) 


for  spirts  oriented  in  the  minus  z  direction,  where  [a.  ,  a.  ]  =  -1.  The  Hamiltonian 
can  be  expanded  as 


+  _ 


+  + 


*  _  - 


K  =  E  +  ELT  a.  a  +iEU  a.  a.  +  |EU  a  a  ..., 
i  i  M  '  i  ij  'i  '  I  ii  ii  ’  i 

®  .  '  * 
where  E  is  the  expectation  value  of  the  energy  in  the  Neel  state,  and  T..  =  T,. 

and  U..  =  U..  .  A  more  symmetric  form  for  is 

i|  |i 


(A3) 


+  - 


+  + 


*  _ 


K  =  E  +ipT..(a.  a.  +  a.  a.  )  +  *EU..a.  a.  +iEU..  a.  a.  -*ET..  . 
0  ij  i|  I  I  I  '  ii  'I  '  I  ij  i|  I  I  i  " 


(A4) 


Thus,  in  this  approximation  the  ground  state  energy  is  different  from  that  in  the  Neel  state 
by  an  amount 


AE®  =  -ilT..  IKo)  1  ,  (A5) 

i  “  n  n 

which  for  a  ferromagnet  is  zero,  of  course. 

In  the  RPA  one  calculates  the  frequency  spectrum  at  finite  temperatures  by  sub¬ 
stituting  for  S  the  thermal  value  of  . .  We  now  show  that  at  low  temperatures  this 
excitation  energy  is  not  exactly  what  one  should  use  In  the  calculation  of  the  free 
energy,  and  that  using  the  correct  frequency  one  no  longer  need  consider  explicitly 
the  zero  point  energy. 

We  start  by  assuming  that  the  temperature  is  low  enough  so  that  we  may  write. 


E  =  E®+AE®+  In.  IKw  .l+iS  n  n  V 

pk  po'  pv  pk /o  k 

*  kk' 


(A6) 


10  35 

This  formula  can  be  mode  to  agree  with  the  first  RPA  or  with  Dyson's  results 

with  the  proper  choices  of  the  coefficients,  V  .  From  this  equation  one 

pk  ,  p'Jc' 

deduces  an  excitation  energy 


Kw.(n)=Kw  +L  Vk'^nk  «'k' 


(AT) 


where  n  stands  for  the  totality  of  occupation  numbers  other  than  n  ,  .  Using  brackets, 
r  pK 

(  )  ,  to  represent  thermal  averages,  one  has 
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where  we  dliHnguish  between  K<i)  l  ^bfch  is  (T  *  0)  and  .  (T)  .  By  expanding 

Pi  pl<  Pi 

the  partition  function. 


2  =  '*■  S  ®  "pk  ^  ,  p'  k'  "pk  "p'  fe'  -  ^ 


(A9) 


one  can  write  the  free  energy  as 


where  F(Km  ,  )  is  the  free  energy  calculated  using  the  zero  temperature  frequencies, 
pk. 

Since 


jiiu 


'pk 


-  »  k  T  — i -  4n(l  -e'  )  *  <n  .  >  , 

®  BKw.l  ^1*. 


pk 


H 


one  has 


o  o  .  ap 

F  =  E  +  AE  +F(K«  ,  )  +  |  E  V.  ..  .<n^.  >^17-; - 

ph.  0P'  Pis. '0  Is.  Pis.  ^^Vk' 

kk' 


(AM) 


(A12) 


or  using  Eq.  A8 

F  -  E°  «  Ae”  *F  (K«.^,)  >  i  E  ,  (Ky  ,,(T)  -  (0)  -Jjf—  ,  (A13) 


dF 


so  that  finally 


F  =  E°  +  PE®+F( 


^%k  ^  ^"pk 


■)  . 


(AM) 


Thus  we  see  that  whereas  is  the  excitation  energy  one  might  observe  in  a 
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resononce  experiment,  the  macroicoplc  propertiei  are  to  be  calculated  using  the  frequency 


36 


These  considerations  are  only  valid  as  long  as  A6  Is  a  good  approx¬ 


imation,  in  which  cose  we  do  not  expect  to  be  able  to  detect  the  thermal  variation  of 
frequency  of  modes  which  are  not  excited,  since  the  zero  point  energy  need  not  be  con¬ 
sidered  explicitly.  In  contrast,  the  zero  point  energy  would  have  to  be  considered  in 
cases  where  the  external  parameten  upon  which  it  depends,  e.g.  the  magnetic  field,  are 
variable. 
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CAPTIONS  FOR  THE  FIGURES 


4. 

5. 


7. 

8. 


The  form  of  the  spin  wove  matrix,  A.. ,  for  YIG.  The  submatrices  are  given  in 
Figs.  3-7. 

The  form  of  the  spin  wave  matrix,  Aj. ,  for  GdIG.  ^  is  the  unit  matrix  and 

420J  .  The  submatrices  are  given  in  Fig.  3-7. 

dc  ac 

The  submatrix  A^^(k) .  Here  Aj  *  ^"^aa  '^"^ac^c' 

Ao  *  -lOJ  cos(kq/4),  and  Ao  ■  -lOJ  co$(3kq/4).  Here  and  in  Fig.  4-9 
i  aa  V  aa 

we  use  the  notation  k  =  k  =  k  =  k. 

X  y  z 

The  submatrix  A^  (k)  .  Here  M|  ®  -5  exp(ikq/8)  and  “  -5  exp(3ikq/8) . 

The  submatrix  Ajj(k).  Here  D  =  -20  +20Jjj  +20Jjj  -  4JgjS^, 

D]  =  "^^dd'  ^2  “  -5Jjjexp{ikq/4),  Dg  =  -5Jjj  exp(ika/2). 

The  submatrix  A  (k).  Here  L,  =  -2{5S  /2)^J  exp(ikoi/8), 
ac  I  c 

k 


oc 


Lj  =  -2(5Sy2)"J^^exp(3ika/8). 

The  submatrix  A^j(k).  Here  G  =  -2{5Sg/2)^J^jexp{ikq/4). 

The  transformed  spin  wave  matrix  of  Eq.  13,  A(k  ,  1).  For  GdIG  the  entire  matrix 
is  solved.  For  YIG  one  solves  the  8x8  matrix  obtained  by  deleting  the  right 
hand  four  columns  and  the  bottom  four  rows.  The  symbols  not  defined  in  the 
captions  of  previous  figures  are  A^  ®  V3  Ag,  Ag  =  2A2  +  Ag,  Mg  =  ^3  Mj  , 

=  Mj  +  Mg,  Lg  *  V3  Lg,  “  L]  +  Lg,  *  -D^  -  Dg,  and  *  -2Dg 

The  transformed  spin  wave  matrix  of  Eq.  13,  A(k ,  ^).  For  GdIG  the  entire  matrix 
is  solved.  For  YIG  one  solves  the  6x6  matrix  obtained  by  deleting  the  right  hand 
four  columns  and  the  bottom  four  rows.  The  symbols  not  defined  in  the  coptions  of 
previous  figures  are  A^  *  ^^^3  ”  •  ^5  ®  +X  Mg,  M^  =  XMg, 

M^  *  XMJ,  Mg  *  M,  +X^m5,  Mg  -  Mg  +X^mJ,  M^  =  Mg  +X^M,  , 
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*■5  ^5  “  ^*■1'  ‘•6  "  ’■2  h‘h  '^^^*■1'  ‘■8  *  •■1  ■^^^'■2' 

L^  ®  L^  +X^L2»  and  D7  =  “^Dg  -  X^D]  . 

10 .  The  spin  wave  spectrum  for  YIG  for  1^  lying  along  a  [111]  direction,  and  for 

J  “  0  and  =  0^  which  corresponds  to  =  20^  cm  '  .  Here 

and  in  the  following  figures  the  polarization  of  the  mode  is  indicated  by  a  plus  or 

minus  sign.  The  eigenvalue  of  the  rotation  operator  is  also  given.  The  modes 

2 

corresponding  to  the  eigenvalues  X  ond  X  are  degenerate. 

1 1 .  The  spin  wave  spectrum  for  YIG  for  k  lying  along  a  [111]  direction,  and  for 

■'ac/-'ad  =  -2  and  =  0,  or  =  26.2  cm"’,  =  5.24  cm  ’,  and 


12.  The  spin  wave  spectrum  for  YiG  for  lying  along  a  [111]  direction,  and  for 

V-'od  "  V-'od  '  ^  ■’=0  *  Jdd  ‘  • 

13.  The  low  frequency  part  of  the  spin  wave  spectrum  of  GdIG  for  1^  lying  along  a 

[111]  direction,  and  for  =  31 .8  cm  \  =  Jjj  =  6.35  cm  ,  and 

J  /j  .  =  .25,  or  J  ,  =  7.00  cm  '  and  J  =  1 .75  cm  V  See  the  text  for 

ac  ac  dc  ac 

a  discussion  of  the  WMF  and  average  frequencies. 

14.  The  low  frequency  part  of  the  spin  wave  spectrum  of  GdIG  for  j(  lying  along 

[111]  direction,  and  for  J  ,  =  31 .8  cm  ,  J  =  J  , .  =  6.35  cm  ,  and 

ad  aa  aa 

J  /J,  =  .0833,  or  Jj  =  3.66  cm  '  and  J  ^  =  0.305  cm  .  For  this  choice 
ac  dc  dc  ac 

of  parameters  it  happened  that  one  mode  was  highly  degenerate  corresponding  to 
eigenvalues  of  the  rotation  operator  indicated  in  the  legend. 

15.  The  magnetization  of  YIG  below  50®  K.  The  dots  are  the  experimental  values  of 

22 

Solt  .  The  curves  were  calculated  for  the  values  of  the  ratios 

Jdd/Jod  indicated  by  the  numbers  in  parentheses;  the  corresponding  values  of 
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16. 


Jod  Jgd  *  17.8  cm"',  (.3,  .1);  «  38j0  cm"’,  and  (.2,  .2): 

=  31 .8  cm  * .  The  effective  magnetic  field  it  given  by  Eq.  37. 

The  magnetic  contribution  to  the  specific  heat  of  GdIG.  The  dots  are  the  ex¬ 
perimental  values^  and  the  curvet  are  calculated  taking  J  ,  =  31 .8  cm  and 
1 

■’aa  *  ■’dd  “  values  of  and  corresponding  to  the  in¬ 

dicated  values  of  J  /J,  are,  J  /j,  *  .25;  J,  *  7.00  cm"’  and 
ac  dc  ac  dc  dc 

J„  “  ^-75  cm  ,  and  J  /J  ■  =  .1667;  J,  ■  474cm  and  J  =0.79  cm 
ac  ac  dc  dc  ac 

The  magnetization  of  GdIG  below  20®  K  taking  =  31 .8  cm  ’, 

Jaa  ~  -^dd  ~  ®”'  '  "^dc  “  Jgg  “  ^  *75  cm  ’ . 


17. 
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